Slow dynamics of a confined supercooled binary mixture II: Q space analysis 
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We report the analysis in the wavevector space of the density correlator of a Lennard Jones 
binary mixture confined in a disordered matrix of soft spheres upon supercooling. In spite of the 
strong confining medium the behavior of the mixture is consistent with the Mode Coupling Theory 
predictions for bulk supercooled liquids. The relaxation times extracted from the fit of the density 
correlator to the stretched exponential function follow a unique power law behavior as a function 
of wavevector and temperature. The von Schweidler scaling properties are valid for an extended 
wavevector range around the peak of the structure factor. The parameters extracted in the present 
work are compared with the bulk values obtained in literature. 
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I. INTRODUCTION 

Modifications of vitrification properties of liquids un- 
der confinement are of interest for relevant technological 
and biophysical problems 1,2 . While experiments show 
that the glass transition is not usually suppressed when 
the glass former is confined, changes in the temperature 
of the transition and in dynamical properties of liquids 
upon supercooling are induced by both the geometrical 
confinement and the interaction with the substrata^. 

In the study of the glass transition phenomenology the 
mode coupling theory (MCT) of the evolution of glassy 
dynamics plays a fundamental role. MCT is able to de- 
scribe the slowing down of supercooled liquids in a re- 
gion of mild supercooling on approaching a temperature 
Tc2i2&. This temperature marks a crossover from a re- 
gion where dynamics is mastered by the transient caging 
of nearest neighbor to a lower temperature region where 
hopping processes dominate. In the ideal version of the 
theory, where hopping is not included, the non linear 
feedback mechanisms in the microscopic dynamics of the 
particles become so strong that they lead to the struc- 
tural arrest of the system. Close to Tc MCT asymptotic 
solutions show analytic behaviors of the density correla- 
tor from which parameters of the theory can be extracted. 
These behaviors have been effectively found in experi- 
ments and simulations on supercooled liquids in regions 
close to Tc where hopping was negligible 5 . 

The role of the cooperative dynamics of the particles 
in the glass transition is still under discussion and there 
is an increasing activity on the definition and observation 
of dynamical heterogeneities^. Cooperativity in the dy- 
namical behavior of the particles increases as the system 
is approaching the glass transition, suggesting the idea 
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of the existence of a dynamic correlation lengthiSaLiLII. 

Molecular dynamics studies of model liquids in re- 
stricted geometries intended to assess the applicabil- 
ity of MCT can give an important contribution to 
the characterization of vitrification processes in confine- 
ment 12 . 13 i 14 . 15 i 16 . 

Liquids confined in network of interconnected pores 
with a large value of porosity, as Silica Xerogels, can 
be appropriately studied with models where the confin- 
ing solid is built as a disordered array of frozen micro- 
spheresiii&iii. We consider here one of such system. A 
liquid Lennard Jones binary mixture (LJBM) composed 
of 80% of particles A and 20% of smaller particles B 
embedded in an off lattice matrix of soft spheres. We 
performed Molecular Dynamics (MD) simulations of the 
confined LJBM upon supercooling in order to test the 
predictions of MCT. 

The bulk phase of this LJBM behaves a la Mode Cou- 
plniifiLsL^. A preliminary analysis on our system 23 
showed evidences that predictions of MCT hold also for 
the confined LJBM. A detailed study on the dynamics 
carried out in the direct space has shown two important 
differences with respect to the bulk upon supercooling 24 . 
The smaller B particles tend to avoid the soft sphere 
interfaces on lowering temperature and correspondingly 
their diffusion coefficient becomes lower than that of the 
A particles, at variance with the bulk where there is 
no such an inversion. Hopping processes are markedly 
present for both A and B particles already above Tc 
while in the bulk system only B particles show consis- 
tent hopping. 

As a possible consequence of these differences, we ob- 
served that the range of validity of the power law fit to 
the diffusion coefficient as predicted by MCT in the long 
time region, late a relaxation region, appears much more 
limited than in the bulk. 

A more refined assessment of the MCT validity can 
be obtained by a quantitative comparison between the 
asymptotic predictions of the theory and the results of 
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our computer simulation in the (Q,t) space, where the 
density correlators of a liquid testing MCT display clear 
signatures of the approach to Tq . In this space important 
scaling relations of the theory can be verified. For these 
reasons we carry on here an analysis of the behavior of 
the self intermediate scattering functions, Fs(Q,t). 

In the next section we recall the predictions of the ide- 
alized version of MCT connected to our analysis. In Sec. 
Ill we describe our model system and Molecular Dynam- 
ics details. In Sec. IV we perform a test of MCT prop- 
erties on our system. In particular we analyze the static 
structure factor and the self intermediate scattering func- 
tion. For this latter quantity we perform numerical fit to 
MCT analytical laws and extract parameters of the the- 
ory. Last section is devoted to conclusions. 



II. MCT SCENARIO FOR THE EVOLUTION OF 
THE STRUCTURAL RELAXATION 



are very close to /g and show a very mild temperature 
dependence 



fQ-fo~ 0(e) 



(2) 



In the asymptotic limit where e << 1 MCT predicts a 
factorization 



t) - /A = h Q G(t) 



(3) 



where \iq is a positive factor which does not depend on 
time. The so called (3 correlator, G(t), obeys the scaling 
relation 



G(t) ~ V\e\9(t/t e ) (4) 
where the time scale t e has been introduced 
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The MCT for the evolution of glassy dynamics devel- 
oped originally to deal with the cage effect provides a 
model for an ideal liquid-glass transition. The evolution 
of the structural relaxation is described in terms of the 
asymptotic behavior of density correlators near singular- 
ities which take place close to a cross-over temperature 
Tc- 

For a density correlator or any other correlator that 
couples with density fluctuations, 4>(Q,t), MCT predicts 
upon supercooling a diversification of relaxation times 
related to the nearest neighbor caging. In the short time 
regime the correlator show a fast relaxation associated 
with the ballistic motion of the particles. The high tem- 
perature correlator relaxes then exponentially for longer 
times. On lowering temperature for intermediate times 
a shoulder appears in the relaxation law. As the tem- 
perature is further lowered the shoulder enhances and 
becomes a plateau and the correlation function clearly 
shows a two step relaxation. The presence of the plateau 
is related to the trapping of the particle in the cage 
formed by its nearest neighbors. The time interval in 
which the correlation function is close-to and above the 
plateau is called the (3 relaxation region. The region that 
starts when the correlation function leaves the plateau is 
called the a relaxation region. 

The long time limit of </>(Q, t) is the non ergodicity 
parameter /q(T). This parameter exhibits a singularity 
at a temperature Tc where the system undergoes an ideal 
transition from an ergodic to a non ergodic behavior. 
/q(T) jumps discontinuously from above Tc to /q at 
Tc- On approaching Tc from below 
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where e = (T — Tc)/Tc is the small parameter of the 
theory and fq is called the critical non ergodicity pa- 
rameter. On approaching Tc from above an effective Jq 
can be identified with the height of the plateau of the cor- 
relator. In this range of temperatures the values of /q 



Here t is a microscopic characteristic time of the sys- 
tem and the exponent a is called the critical exponent. 
In the time scale which corresponds to the decay of the 
correlator towards the plateau, t <t c 

g(t/U) = {t e /t) a (6) 

For the correlator decay of the liquid from the plateau 
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which is valid for times larger than the time scale of the (3 
relaxation and much smaller than the a-relaxation time 
t, t e < t << t. The exponent b is known as the von 
Schweidler exponent. The critical and the von Schweidler 
exponents are related by the formula 



A = 



[r(i- a )] 2 _ [r(i + b)f 



r(l-2a) r(l + 26) 



(8) 



where A is called the exponent parameter and r(x) is the 
Euler Gamma function. 

In the late part of the a-relaxation regime the shape 
of the correlator of a liquid approaching the MCT 
crossover temperature^ can be represented by means of 
the Kohlrausch- Williams- Watts (KWW) analytical func- 
tion 



<t>(Q,t) = f Q exp(-(t/Tf) 



(9) 



where /3 is the Kohlrausch exponent. 

The a-relaxation time r is expected from MCT to fol- 
low the power law 



T- L ~(T-T c y 



(10) 



Both the exponent 7 and the crossover temperature Tc 
should coincide with the ones appearing in the tem- 
perature dependence of the diffusion coefficient D ~ 
(T-TcV- 



An important prediction of MCT is the validity of the 
time-temperature superposition principle. It states that 
in the low temperature regime, approaching from above 
the MCT crossover temperature, it is possible to time- 
rescale the correlators evaluated at different T into a sin- 
gle master curve: 

m = 4>(t/ T (T)) (ii) 

4> is the master function which in the late a-relaxation 
limit becomes the KWW function 

In the late /3, early a relaxation regime Eq. [3] can be 
written as 

HQ,t) = f c Q -h Q (t/ T ) b (12) 

This equation is called the von Schweidler law and r in- 
dicates a characteristic time of the system. 

The exponent 7 in Eq. is related to the exponents 
a and b by the equation 



While both t and the exponent (3 are Q dependent, 
the exponent b should not depend on Q. These two ex- 
ponents are related: j3 > b and f3(q — > 00) = 




FIG. 1: Static structure factors. The temperatures corre- 
sponding to different curves are 5.0, 2.0, 0.8, 0.58, 0.538, 0.48, 
0.465, 0.43, 0.41, 0.39, 0.37. For clarity we shifted up the lower 
temperature plots. For AA the shift is 0.2, for BB is 0.025 
and for AB is 0.1. 



III. MODEL AND MOLECULAR DYNAMICS 

The model studied consists in a rigid off lattice matrix 
of 16 soft spheres where a liquid ^80^20 Leonard Jones 
binary mixture (LJBM)^!* 2 ^ is embedded. The mixture 
is composed of 1000 particles. In the following LJ units 
will be used. The parameters of the LJBM are caa = 1, 
oaa = 1, ess = 0.5, obb = 0.88, 6ab = 1.5, and 
gab — 0.8. The parameters of the soft sphere poten- 
tial are esA = 0.32, 0-5,4 = 3, e$B = 0.22, <jsb = 2.94. 
The box length is L=12.6. We performed Molecular dy- 
namics simulations in the NVE ensemble. Wc studied a 
set of temperatures ranging from T=5.0 down to T=0.37. 
The timestep for the lowest temperature is 0.02 and the 
largest production run consists in 14xl0 6 timesteps. For 
all temperatures the time range investigated by MD was 
sufficient for the correlation functions to decay to zero in 
the long time limit. Further MD simulation details are 
reported in ref.sS&Si. We did not find any qualitative dif- 
ferences in the behavior of the thermodynamic quantities 
with respect to the bulk upon supercooling^*. 

IV. ANALYSIS OF THE DYNAMICAL 
BEHAVIOR IN (Q,t) SPACE 

A. Static structure factor 

In the framework of MCT the static structure factor 
is the input information for MCT equations^. This 




FIG. 2: Intermediate scattering function for A particles (up- 
per panel) and for B particles (lower panel) for all the in- 
vestigated temperatures. The continuous lines are calculated 
from the MD trajectories at the peaks of the structure factor 
Qmax — 7.0 and Qmax — 5-8. The long dashed lines are 
the fit to the KWW stretched exponential function, see Eq|5| 
In the lower graph the dotted lines are the A particles data, 
reported for comparison. 



function is expected to depend smoothly on wavevectors 
and control parameters and it is also expected to show 
no singularities upon supercooling. In Fig. ^ we report 
the structure factors S(Q) of the confined mixture for 
the temperatures investigated. They behave consistently 
with MCT as the functions smoothly vary upon super- 
cooling implying no divergence of correlation length. No 
signature of phase separation of the mixture are evident 
in confinement as we already deduced from the analy- 
sis of the gAB{ r ) 2i - We observe that there are no rel- 
evant discrepancies in the confined SV/(Q) with respect 
to the bulkSS neither in the shape nor in the peak po- 
sitions. For Q values around the first structure factor 
peak MCT features are best evident in the density cor- 
relators. The peak positions in our confined LJBM show 
a very mild temperature dependence, especially for low 
temperatures. We will therefore consider them constant. 
The values Q^ax ~ 

7.06 

Qmax = 5-90 will be used to 
test MCT as a function of temperature in the following. 
In the case of the bulk the values used were respectively 
Qmax = ^-^5 and Q^ax = 5-75. 



B. Self intermediate scattering function 

We consider now the self intermediate scattering func- 
tion (SISF) defined as 

Fs (Q ; t ) = -<J2 e-^^We^W > (14) 




FIG. 3: In the main frame, values of the relaxation time r 
obtained from the fit to the KWW of the SISF (see Fig|5J 
for A (empty circles) and B (empty squares) particles. In 
the same graph values of the diffusion coefficients D obtained 
from the MSD2£ for A (+) and B (X) particles. The power 
law fits are also plotted. These fits are performed on the 
temperature window 0.41 < T < 0.58, limited by vertical 
dotted lines. Dashed lines are the fit for A and continuous 
lines are the fit for B particles. Parameters obtained from the 
fit are reported in Tab[I] In the inset /q for A (filled circles) 
and B (filled squares) particles and stretching exponent j3 
for A (empty circles) and B (empty squares) particles are 
reported. 



This correlator can be used to test the predictions of 
MCT close to T c as described in Sec. II. In Fig. □ the 
SISF of our confined system is plotted for all considered 
temperatures and for both particle types at the peak of 
the structure factors. At high temperatures the corre- 
lators relax exponentially for longer times. On lowering 
temperature we observe the onset of the plateau related 
to the two step relaxation typical of a supercooled liquid. 

An onset temperature for the plateau of T w 0.6 can 
be identified for our system. The corresponding temper- 
ature for the bulk is T « 1.0. Analogous to bulk super- 
cooled liquids for the lowest temperature we observe that 
correlators stretch several decades in time. We also see 
that A and B particles relax in a similar fashion but larger 
A particles relax faster at this wavevectors, especially fol- 
lower temperatures. This counterintuitive phenomenon 
that happens only in confinement is possibly connected 
to packing constraints. We have in fact detected that on 
lowering temperature smaller B particles tend to avoid 
the soft sphere interfaces and preferentially move in the 
center of the interstices. Correspondingly diffusion co- 
efficients of A and B particles extracted from the mean 
square displacement (MSD) show an inversion as super- 
cooling progressesSi. 
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FIG. 4: Wave vectors analysis at T = 0.41. The range for 
both particles is 2.5 < Q < 14.5, with a step of 5Q — 1.5. 
The KWW fit (long dashed lines) is also shown. 
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FIG. 5: In the graphs are shown the KWW parameters fg, r, 
f3 upon varying the wave vector Q, for A (panels on the left) 
and B particles (panels on the right). The various symbols 
refer to different temperatures: circle 0.37, square 0.39, dia- 
mond 0.41, triangle up 0.43, triangle left 0.48, triangle down 
0.58. The dotted line of the bottom panels refers to the val- 
ues extracted from the von Schweidler fit 6a = 0.355 and 
b B = 0.350. 



C. Kohlrausch- Williams- Watts analysis of the self 
intermediate scattering function 

In Fig. together with the SISF we show the fit to 
the KWW function, as defined in Eq. performed for 
temperatures T < 0.8 in the a relaxation region for both 
A and B particles. A remarkable agreement with the data 
is found. From the fit the KWW relaxation time t, the 
stretching exponent (3 and the effective non ergodicity 
parameter Jq are extracted. 

MCT predicts both D and t" 1 to follow the same 
power law, see Eg 1101 The relaxation times r are shown 
in the main frame of Fig- EI In the same picture we show 
the diffusion coefficients D previously reported and ob- 
tained from the slope of the MSD24. We also display in 
the figure the fit to the MCT predicted behavior. Values 
of Tc and 7 obtained from the fit are reported in Tab[3 
We observe that both in bulk and confined system the 
crossover temperature Tc extracted from D and r and 
from A and B particles coincide, Tc — 0.36 A substan- 
tial reduction of Tc is observed in confinement. The bulk 
value reported in literature is Tc — 0.43. 

The exponents 7 which determine the slope of the fit 
are the same for A and B particles both for r and D. This 
agrees with the prediction of MCT that these exponents 
should be universal and do not depend on the species 



of the particles involved. A discrepancy with respect to 
MCT prediction is represented by the different value of 7 
between D and r. As we can see from the table the same 
discrepancy affects the mixture also in the bulk phase and 
it may be due to long time activated processes causing 
a breakdown of the predictions of the theory only for 
quantities calculated at extremely long tiniest. In fact 
in confinement, where hopping is consistent for both A 
and B particles, a reduction of circa 35% is observed for 
7 in going from t to D. In the bulk phase, where hopping 
is more consistent for B particles, the reduction observed 
is 30% for B and 20% for A particles. 

A substantial shrinkage of the range of validity of the 
theory is found in confinement. In Fig. [21 we in fact see 
that the MCT predicted power law behavior is verified 
for t in the temperature window 0.58 < T < 0.41, the 
same found for the D coefficients^, which corresponds 
to 0.153 < e < 0.631 while in the bulk 0.07 < e < 1.30 20 . 

The values of the parameters /q and (3 extracted from 
the fit are plotted as a function of temperature in the 
inset of Fig. [3J We note that the effective non ergodic- 
ity parameter /q only slightly depends on temperature 
in agreement with MCT that predicts a consistent tem- 
perature dependence of this factor only below Tc- The 
stretching parameter (3 decreases upon lowering the tem- 
perature as also expected, but we observe no flattening of 
this exponent on approaching Tc as found for bulk sys- 
tems. Besides, correlators appear much more stretched in 
confinement. In the bulk (3 > 0.78 while here the lowest 
values are slightly below 0.4. 

We report in Fig. 0] the temporal evolution of the cor- 
relators for different wavevectors Q ranging from 2.5 to 
14.5 for the temperature T = 0.41. In the same picture 
we also show the fit to the KWW function. A remarkable 
agreement is found. The height of the plateau strongly 
depends on Q. The same analysis was performed for 
T = 0.58, 0.48, 0.43, 0.39 (not shown). 

The values of the parameters of the fit are reported for 
all the above temperatures in Fig. All the parameters 
show a monotonic behavior in Q. In particular at high 
temperature the relaxation time r spans two decades. It 
is worth noting that the f3 exponent is converging, for 
T — ► Tc and for Q — > 00 to the value of b for both 
types of particles as expected^. Finally we note that 
the prefactor Jq is slightly T dependent for all Q. 

In Fig. [HI the r values extracted from the fit are re- 
ported as a function of the shifted temperature T — Te- 
la the temperature range where the power law behavior 
is effective, that for our system is 0.41 < T < 0.58, 7 is 
expected to be independent on Q. In Fig.EJwe effectively 
observe that the curves are parallel to each other and par- 
allel to the straight line overlayed that corresponds to the 
value of 7 = 2.9 extracted from the fit performed at the 
peak of the structure factor (see Tab^). 




FIG. 6: Relaxation times r (symbols) vs. shifted temperature 
T — Tc- The power law fit for Q = Qmax is also plotted. Var- 
io 
st 




t/T 

FIG. 7: Intermediate scattering functions of Fig|5] rescaled 
with respect to the corresponding KWW relaxation time for 
both A (top panel) and B (bottom panel) particles. The 
dashed curves represent the von Schweidler fit. 



D. von Schweidler test for the self intermediate 
scattering function 

According to Eq. ^2 the time dependence of the cor- 
relation functions can be rescaled by the relaxation time 
r extracted from the fit to the KWW. The curves for 
the lowest temperatures are expected to collapse onto a 



FIG. 8: Values of the parameters of the von Schweidler fit 
(black diamond, white circles) and for comparison the pa- 
rameters of the KWW (black triangles, white squares) upon 
varying Q for A particles. The static structure factor is rep- 
resented by the dotted line at T = 0.37. 

mastercurve in the late /3, early a relaxation region. In 
Fig. |7| we report for both particle types the master plot 
of the SISF of Fig. In the same figure we show the fit 
to the von Schweidler functional form of Eq. El 

The values of the parameters obtained for A and B 
particles from the fit are: 

Iq.a = 0-72, h Q , A = 0.35, b A = 0.355 

Jq.b = 0-79, h Q , B = 0.44, b B = 0.350 

In the bulk the values of the exponent b are respectively 
0.51 for A and 0.46 for B particles while we observe that 
in confinement b values are lower and almost the same. 
This might indicate that hopping is influencing also the 
earlier a relaxation region. 

An important prediction of MCT is that the von 
Schweidler b value should be independent of the wave 
vector Q. This test was performed in the bulk22, where 
the authors stressed a small dependence of b between 
Qmax and Qmin, where Qmin is the position of the 
first minimum after the peak of the static structure fac- 
tor. In Fig. El the values for /q and b are presented as a 
function of Q for A particles. The values for /q and 
for the KWW and the S(Q) are also shown. The Jq fac- 
tors are in perfect agreement with the Gaussian behavior 
predicted by the MCT and found in the bulk simula- 
tions. We can therefore consider the /q KWW factors as 
a good approximation of the critical non-ergodicity pa- 
rameter, /q, at such low temperature. The b parameter 
decreases monotonically, except between the first S(Q) 
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maximum and the following minimum where it is practi- 
cally constant, similar to the bulk. We note that the j3 
values are always higher than the relative b values, this 
is a consequence of the fact that the KWW can describe 
the late a but not the late /?, early a regime where the 
data are well fitted by the von Schweidler law. 

In our correlators oscillations, possibly due to finite 
size effects, in the time zone where the power law of Eq|I)] 
holds, prevents us from calculating the exponent a of 
MCT directly. It is still possible to evaluate the param- 
eter a from the closure relation of Eq. ^| Using the 
7 calculated from r and the von Schweidler exponent b 
found previously, the values of the critical exponent a for 
the two species are respectively 

a A = 0.335, a B = 0.342 

These values are quite high, although within the topmost 
value indicated by the theory a < 0.395. 

V. SUMMARY AND CONCLUSIONS 

We have studied the dynamics of a model glass former 
in confinement. This study is intended to assess the valid- 
ity of MCT for supercooled liquids embedded in strongly 
repulsive and highly confining environments. Our model 
consists in a LJBM embedded in a disordered matrix of 
soft spheres. It represents a possible modellization for 
liquids hosted in Silica Xerogels. 

The analysis performed in r space and previously re- 
ported^ showed that the MSD behavior can be inter- 
preted for both A and B particles in terms of the onset 
of the cage effect and that the diffusion coefficients D ex- 
tracted from the MSD can be fitted with the power law 
predicted by MCT. There is however a clear evidence that 
the range of validity of the MCT predictions in the late a 
relaxation region suffers a reduction of 60% with respect 
to the bulk. In terms of the small parameter of the theory 
e the bulk range is 0.07 < e < 1.30 while in the confined 
system it results to be 0.153 < e < 0.631. We inferred 
that the upper bound reduction might be due to the fact 
that B particles tend to avoid the soft sphere interfaces 
in lowering temperature. The MCT power law starts in 
fact to hold for D only when most B particles are avoid- 
ing the soft spheres^. The lower bound reduction is to 
be connected to the presence of important hopping pro- 
cesses more marked with respect to the bulk, which tend 
to introduce deviations from the ideal MCT predictions 
earlier than in bulk. 

We have presented here an analysis in (Q,t) space of 
the confined LJBM by considering the SISF and we have 
tested the main scaling relations of MCT. 

The static structure factor is similar to that of the 
bulk phase. The SISF is well fitted with the stretched 
exponential function and the power law behavior of the 



a-relaxation times r extracted for both A and B particles 
is consistent with MCT. Power law fit to both D and r 
as a function of temperature give the same value for the 
crossover temperature Tc — 0.355. We observe a reduc- 
tion of circa 20% of Tq in going from bulk to confined 
which can be ascribed to the existence of an upper bound 
in confinement for the size of the domains of cooperative 
dynamics. A discrepancy observed in the values of the 
exponents 7 obtained for D and r is to be connected 
to the more important presence of hopping processes in 
confinement. Hopping processes influence in fact the be- 
havior of D more than that of t causing a decrease of the 
value of the exponent. This discrepancy was present also 
in the bulk. 

The existence of a confining disordered host structure 
seems to be connected also to the lowering of the stretch- 
ing parameter possibly caused by a larger distribution of 
relaxation times with respect to the bulk. 

The behavior of the SISF as a function of Q is very sim- 
ilar to that of the bulk. All parameters show a monotonic 
behavior as a function of Q and the stretching exponent 
approaches b for high Q. The exponent 7 shows no Q 
dependence. 

The SISF plotted against t/r shows a scaling behavior. 
The mastercurve can be fitted with the von Schweidler 
function in the late j3, early a region. Analogous to the 
bulk, the exponent b is practically constant between the 
first maximum and the following minimum of the struc- 
ture factor. Similar to /3 a substantial reduction of b is 
observed in going from bulk to confined LJBM. 

The exponent parameter A of the system can be cal- 
culated from Eq. |S| Contrary to MCT prediction we get 
different values of A on using the exponent a or the ex- 
ponent b in Eq. [S] 

\ a , A = 0.680, A M = 0.871 

\ a ,B = 0.663, X b:B = 0.873 

This inconsistency is present also in the bulk but it is 
much less severe. The value of A reported for the bulk is 
0.79 ±0.04. Here we have a much larger uncertainty that 
leads to 0.8 ±0.1. 

Considering the strong confinement experienced by the 
LJBM, the agreement between MCT predictions and the 
behavior of our confined system as a function of T and Q 
is most remarkable. Differences in behavior with respect 
to the bulk on approaching Tc are in fact to be ascribed 
to the presence of more relevant hopping effects. 

The possibility to connect the decrease of Tc and the 
enhancement of hopping effects in confinement with a 
key role played by dynamical heterogeneities in the glass 
transition scenario is an issue that would deserve to be 
addressed in future studies. 
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TABLE I: Values of the parameters of Eg . 1101 extracted from 
the fit to r and D. Values for the bulk system are from 
RefsM. 





mA rrB A B 
T c T c 7 7 


Confined system 




From t 
From D 


0.355 0.355 2.90 2.89 
0.356 0.356 1.86 1.89 


Bulk system 




From t 
From D 


0.432 0.432 2.6 2.6 
0.435 0.435 2.0 1.7 
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